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RESTORATION  OF  ATMOSPHERICALLY  DISTORTED  IMAGES 


PROGRESS  REPORT 
by 


James  L.  Harris 


1.0  STATEMENT  OF  THE  PROBLEM 


The  space  age  has  created  renewed  interest  in  the  problems 
associated  with  optical  observations  from  the  ground  of  objects  in 
space.  It  is  possible  to  build  large  optical  systems  whose  image 
quality  is  limited  only  by  diffraction.  In  practice  such  ideal 
imagery  is  not  obtained  because  of  the  image  deterioration  induced 
by  atmospheric  turbulence. 

The  Visibility  Laboratory  has  been  conducting  a  program  of 
research  on  possible  means  of  accomplishing  restoration  of  these 
atmospherically  distorted  images.  This  work  has  been  sponsored 
by  the  Advanced  Research  Projects  Agency  by  transfer  of  funds  to 
this  Laboratory^  existing  Bureau  of  Ships  contract  NObs-84075. 
Considerable  progress  has  been  made  in  the  development  of 
techniques  which  may  be  suitable  for  the  restoration  of  these 
images.  This  report  summarizes  the  effort  on  this  project  from 
its  conception  in  December  of  I960  through  December  1962. 

The  contents  of  this  report  have  been  placed  in  more  or  less 
historical  order  because  the  present  day  concepts  have  evolved 
from  the  earlier  work  and  it  is  felt  that  the  presentation  in  this 
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order  of  time  may  make  it  e&eier  for  the  reader  to  follow  the 
development. 


2.0  SIMPLE  MATHEMATICAL  DISTORTIONS  AND  RESTORATIONS 

2.1  Introduction 

At  the  time  that  the  project  was  initiated,  specific  image 
restoration  techniques  suitable  for  this  application  had  yet  to 
be  developed.  It  was  felt,  however,  that  advances  in  technology 
in  fields  such  as  information  theory,  communication  theory,  and 
related  sciences  made  it  a  reasonable  time  to  take  a  new  look  at 
this  very  old  problem. 

Atmospheric  turbulence  can  be  said  to  be  a  form  of  noise 
acting  on  the  image.  Processes  of  extraction  of  signals  in  the 
presence  of  noise  have  long  been  used  in  radars  and  other  sensors. 
Techniques  of  integrating  many  pulses  in  order  to  improve  the 
signal-to-noise  ratio  have  been  utilized  in  radar  systems  since 
their  development  during  World  War  II.  This  form  of  extraction 
of  signal  from  noise,  i.e.,  integration  of  signal  and  noise,  is 
applicable  in  the  radar  case  because  the  noise  is  additive. 

That  is,  the  noise  serves  to  amplitude  modulate  the  signal.  In 
the  case  of  the  image  distortions  due  to  atmospheric  turbulence, 
the  noise  serves  to  smear  the  image  spatially  and  it  would 
therefore  be  expected  that  the  direct  integration  techniques 
applicable  to  amplitude  modulated  noise  conditions  must  be 
modified  for  this  application. 
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2.2  An  Example  of  a  Simple  Distortion  and  Restoration 

A  very  early  mathematical  exercise  of  familiarization  with 
spatial  distortions  may  serve  to  demonstrate  one  way  in  which  a 
series  of  distortions  may  be  utilized  to  extract  information  by 
which  a  single  image  can  be  constructed  whose  image  quality  is 
superior  to  that  of  any  of  the  individual  distortions.  Figure  1 
shows  a  sample  block  letter  and  one  possible  distortion  of  the 
letter.  The  distortion  is  obtained  by  allowing  each  square 
forming  the  block  to  move  left,  right,  up  or  down,  each  move 
with  equal  probability.  In  the  case  of  this  simple  model,  it 
was  assumed  that  each  square  always  moves,  that  is,  the 
probability  that  a  square  will  remain  in  its  original  position  is 
zero.  This  simple  model  assumed  that  there  were  only  two  levels 
of  intensity  possible .  If  two  elements  from  the  original  block 
letter  moved  into  the  same  square,  that  square  would  become  black 
and  would  not  be  distinguishable  in  level  from  a  square  which 
contained  only  one  element  from  the  original  block  letter.  The 
manner  in  which  the  distorted  image  is  formed  may  be  observed  by 
noting  the  arrows  in  Fig.  1  showing  the  movement  of  each  of  the 
original  elements  of  the  block  letter.  Figure  2  shows  17  dis¬ 
tortions  of  an  arbitrarily  selected  block  letter.  Based  on 
observations  of  the  17  distortions,  it  is  not  apparent  which 
block  letter  is  actually  present. 

For  this  type  of  distortion,  restoration  can  be  accomplished 
by  what  has  been  termed  nforbidden  image  diagrams.”  If,  in  the 
distortion,  we  can  find  a  blank  square  which  is  bounded  above, 
below,  to  the  left  and  to  the  right  by  blank  squares,  then  we  can 


Figure  I.  Sample  block  letter  and  a  simple 
spatial  distortion 


SIO  Ref.  63-10 


■ 

■ 

■ 

■ 

■ 

■ 

s 

H 

HI 

mmm 

Ei 

I 

mmm 

_ 

|9B 

mmm 

MB 

w/m/M 

SM 

tmwM, 

c 

w Bp 

4w/mM 

Ip 

I  2  3  4  5  6  7  8 


_ _ _ 

1 

rift 

It  ^ 

|f^Jp 

■ 

j§ 

Ip 

s 

9 

10 

1  1 

12 

13 

14 

15 

16 

Figure  2.  A  simple  spatial  distortion  and  recovery 
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deduce  that  the  original  undistorted  image  was  blank  in  the  central 
square.  This  follows  directly  from  the  assumptions  of  the  distortion 
since  if  the  central  square  had  contained  an  element!  it  is  forced 
to  more  up,  down,  left  or  right*  and  hence  would  show  up  in  the 
distortion  in  one  of  these  squares.  The  "forbidden  image  diagram" 
for  a  single  distortion  is  therefore  a  map  which  shows  those 
increments  which  could  not  hare  contained  an  element  of  the  target, 

A  "forbidden  image  diagram"  can  be  made  for  each  of  the  17 
distortions.  Restorations  are  accomplished  by  noting  that  the 
"forbidden  image  diagrams"  are  cumulative;  that  is*  we  can  super¬ 
impose  each  of  the  "forbidden  image  diagrams"  to  obtain  the  sum 
total  of  our  information  about  those  squares  which  could  not  have 
contained  a  target  increment.  Figure  2  also  shows  the  construction 
of  a  cumulative  "forbidden  image  diagram,"  The  "forbidden  image 
diagram"  shown  below  distortion  2  contains  those  squares  for  which 
a  target  element  could  not  be  present  as  deduced  from  both 
distortion  1  and  distortion  2,  The  "forbidden  image  diagram"  below 
distortion  17  therefore  contains  all  squares  which  are  known  to 
have  not  contained  a  target  element  as  deduced  from  all  17  distortions. 
As  may  be  seen  from  the  figure,  the  block  letter  S,  which  was  the 
block  letter  used  in  the  distortions,  has  been  successfully  recovered. 

The  example  which  has  been  described  was  educational  in  that  it 
demonstrated  that,  at  least  for  this  simple  type  distortion,  there 
were  techniques  by  which  spatial  distortions  could  utilise  informa¬ 
tion  from  a  large  number  of  distortions  to  create  a  single  image  of 
quality  superior  to  that  of  any  of  the  individual  distortions.  The 
model,  however,  is  not  a  realistic  one  in  terms  of  the  type  of 
distortion  to  be  found  as  the  result  of  atmospheric  turbulence. 


3*0  INTRODUCTION  TO  SPATIAL  PHEQUENCIES 


The  later  techniques  of  restoration  employ  the  concepts 
of  spatial  frequencies  and  an  introduction  to  these  oonoepts 
will  be  made  at  this  time* 

Any  function  having  only  a  finite  number  of  discontinui¬ 
ties  can  be  expressed  over  a  defined  region  by  an  infinite 
series  of  sine  and  cosine  terms*  Thus  the  Fourier  expansion 
of  f(x)  over  an  interval  X  is 

f(x)  =  A0  +  A,cos2rTfx  +  A2Cos4Ttfx  +  AjCos6iTfx  +  ... 

r 

+  B,  sin2rrfx  +  B2Bin4jrfx  +  B^sin6tTfx  +  • . . 


where  A^,  B^,  B2,  ...  are  the  sine  and  cosine 

coefficients  which  are  defined  by  Eqs*  (2)  and  (3). 


f (x) cos(i2nfx)dx 


Bi  = 


+x 

P2 


4 


f  (x)  sin(i2rTfx)dx 


This  type  of  analysis  has  become  widely  known  through  its 
application  to  the  analysis  of  electrical  circuit  problems, 
where  a  signal  is  described  in  terms  of  its  frequency 
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composition,  The  Importance  of  these  oonoepta  has  been  amply 
demonstrated  in  the  case  of  radar,  communications  systems, 
and  audio-amplifier  design.  In  practical  problems  an  infinite 
series  is  seldom  involved  and  it  is  generally  recognised  that 
an  approximation  to  a  function  can  be  made  by  including  a  finite 
number  of  terms  with  the  degree  of  approximation  dependent  upon 
the  number  of  such  terms  which  axe  employed.  There  are  obviously 
no  restrictions  on  the  dimensions  of  the  function  which  is  to 
be  expanded  in  Fourier  series.  For  example,  the  function  might 
be  voltage  as  a  function  of  time,  average  annual  rainfall  as  a 
function  of  latitude,  or,  as  in  the  case  of  interest  here, 
image  intensity  as  a  function  of  a  spatial  dimension. 

Consider  the  case  of  a  piece  of  film  on  which  an  image  is 
stored,  A  complete  description  of  the  image  involves  the 
defining  of  the  transmission  of  the  film  as  a  function  of  the 
X-Y  coordinates  of  the  film.  If  X  is  a  horizontal  dimension  then 
it  can  be  imagined  that  a  trace  in  the  X  direction  for  a  fixed  I 
will  result  in  a  function  of  transmission  as  a  funotion  of  X. 
Since  it  is  necessary  to  represent  the  image  only  over  the  film 
format,  this  transmission  as  a  function  of  X  can  be  expanied  in 
a  Fourier  series.  The  coefficients  for  this  expansion  will  be 
dependent  upon  the  value  of  I  which  was  selected  for  the 
trace,  that  is,  each  of  the  coefficients  is  a  function  of  I* 

The  coefficient  itself  can  therefore  be  expanded  in  a  Fourier 
series  as  a  funotion  of  Y,  If  this  is  accomplished,  as  has  been 
done  in  Appendix  I,  the  result  is  a  two-dimensional  Fourier 
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eeriee  which  describee  the  transmission  of  the  film  as  a  function 


of  X  and  I  in  terms  of  an  infinite  series  of  sine  and  cosine 


terms.  One  form  of  this  two-dimensional  Fourier  series  is 


oo  oo 

f(x,y)  =  A  +  2  2  A.  .oos(2nif  x  +  2rrjf  ▼) 

0  i=0  j=0  1J  x  r 


OO  00 

+  2  2  B,,sin(2TTifxx  +  2njfvy)  (4) 

i=0  »  1J  * 


Each  of  the  terms  in  this  series  expansion  is  a  spatial 
sinusoid.  Figure  3  is  a  sketch  of  one  such  sinusoidal 
component.  The  terms  i  and  j  indicate  the  number  of  cycles  which 
the  sine  wave  makes  in  the  X  and  Y  directions  respectively; 
thus  in  the  example  shown,  i  =  3  and  i  =  4»  One  reason  that 
the  use  of  Fourier  expansion  has  particular  significance  in 
optical  systems  is  that  an  optical  system  passes  no  spatial 
frequencies  above  a  defined  cutoff  frequency,  A  finite  series 
is  therefore  all  that  is  required  to  completely  define  a 
bounded  optical  image.  It  is  shown  in  Appendix  II  that  the 
statement  of  a  cutoff  frequency  for  the  optical  system  is 
the  Fourier  equivalent  of  the  usual  statement  with  reference 
to  the  angular  resolution  associated  with  a  diffraction 
limited  optical  system.  Because  of  this  spatial  frequency 
cutoff  in  any  real  optical  system,  it  is  possible  to  define 
any  real  image  by  a  finite  set  of  numbers.  The  significance 
of  this  fact  may  be  better  appreciated  after  the  develop¬ 
ment  of  those  restoration  techniques  which  employ  the  concepts 
of  spatial  frequencies. 
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4.0  A  GEOMETRIC  DISTORTION  MODEL 

4.1  Description  of  the  Model 

If  a  distant  star  is  observed  with  an  optioal  system  having 
a  small  entrance  pupil  diameter,  the  primary  effect  of  the 
atmospheric  turbulence  is  to  move  the  image  about  in  a  random 
fashion.  This  can  be  explained  on  the  basis  of  the  turbulence 
creating  index  of  refraction  gradients  before  the  entrance 
pupil  of  the  optical  system,  thus  acting  in  the  manner  of  a 
prism.  These  simple  observations  serve  as  a  basis  for  a 
geometric  model  of  image  distortion  due  to  atmospheric  turbu¬ 
lence.  In  this  model  it  is  assumed  that  the  entrance  pupil  can 
be  quantized  into  sectors,  each  sector  having  a  refractive 
index  gradient  which  is  determined  in  some  statistical  fashion. 
Figure  4  shows  a  cross  section  through  an  optical  system, 
showing  the  manner  in  which  the  different  sectors  of  the  lens 
form  images  of  the  object  displaced  with  respect  to  one  another. 
The  addition  of  these  independent  images  is  the  resultant 
distorted  image  as  predicted  by  this  geometric  model. 

In  order  to  more  fully  understand  the  nature  of  the  dis¬ 
tortion  process,  a  single  spatial  frequency  will  be  distorted 
in  the  same  manner  as  indicated  in  Fig.  4*  Figure  5  shows  the 
separate  images  formed  by  each  sector  of  the  entrance  pupil  of 
the  optioal  system  and  again  schematically  shows  that  each 
sector  forme  an  image  of  the  particular  spatial  frequency  with 
the  index  of  refraction  gradient  causing  the  image  from  each 
sector  to  be  displaced  right  or  left  from  the  correct  image 
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Figure  4 


Schematic  illustration  of  the  nature 
of  the  geometric  distortion. 


Figure  5  Geometric  distortion  of  a  single 
spatial  frequency. 
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position*  The  composite  image  is,  as  before,  the  sum  of  the 
images  from  each  section  and  in  the  case  of  the  addition  of 
the  sinusoidal  components,  the  resultant  is  also  a  sinusoid 
whose  amplitude  will  always  be  less  than  the  amplitude  of  an 
undistorted  image,  that  is,  the  addition  of  sinusoids  of 
different  phases  results  in  a  sinusoid  haying  an  amplitude  less 
than  would  be  obtained  had  the  phases  been  equal  and  haying  a 
phase  shift  with  respect  to  the  phase  of  the  undistorted 
image.  Thus,  the  geometric  model  predicts  that  in  terms  of 
spatial  frequencies  the  effect  of  atmospheric  turbulence  is 
to  attenuate  the  amplitude  of  the  spatial  frequency  and  shift 
its  phase. 


4.2  An  Example  of  a  Geometric  Distortion 

The  distortion  model  described  in  the  preceding  paragraph 
was  used  to  study  possible  restoration  techniques.  A  necessary 
first  step  in  the  application  of  the  model  is  the  assumption  of 
a  probability  distribution  for  displacement.  As  a  matter  of 
mathematical  convenience,  the  displacement  distribution  or 
gradient  distribution  was  chosen  to  be  a  triangular  distribu¬ 
tion  as  shown  in  Fig.  6.  Sinoe  the  distribution  has 
symmetry  about  zero  gradient,  the  average  displacement  is  zero. 

Figure  7  shows  the  undistorted  image  used  in  the  example. 
Mathematically  the  image  is  &  pulse  as  represented  by  the  first 
ten  Fourier  coefficients.  By  analogy,  it  might  be  imagined 
that  Fig.  7  represents  a  horizontal  trace  through  the  Image  of 
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Figure  6  Probability  distribution  of  the 
refractive  gradient. 
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IMAGES  AND  THEIR  SPECTRA 


Figures  7  through  14 


SIO  Ref.  63-10 


-17- 


a  vertical  bar.  The  representation  by  a  limited  number  of  Fourier 
coefficients  is  &  rough  analogy  to  the  diffraction  limit  because, 
as  stated  in  Section  3.0,  a  diffraction  limited  optical  syetem 
passes  no  frequencies  above  a  definable  cutoff  value. 

The  image  of  Fig.  7  was  subjected  to  the  distortion  process 
illustrated  in  Fig.  4  using  the  gradient  distribution  of  Fig.  6. 
The  lens  was  assumed  to  be  divided  into  ten  segments  and  by 
obtaining  ten  samples  from  the  gradient  distribution,  sample 
distortions  were  constructed.  .  Figures  8  and  9  are  examples 
selected  from  a  group  of  30  such  distortions. 

It  was  noted  in  Section  4«1  that  the  effect  of  this  geometric 
distortion  is  to  attenuate  and  shift  the  phase  of  any  given 
spatial  frequency.  This  observation  leads  to  one  possible 
restoration  technique  which  will  now  be  illustrated.  Figure  10 
shows  the  spectrum  of  the  undistorted  image.  The  dashed  line 
is  the  envelope  and  the  vertical  lines  are  the  discrete 
frequencies  of  the  Fourier  series  which  represent  the  function 
over  the  extent  of  the  image.  Figures  11  and  12  show  the  spectra 
of  the  two  sample  distortions  and  indicate  that  the  amplitude  of 
any  given  spatial  frequency  varies  from  distortion  to  distortion. 
Since  it  has  been  observed  that  a  spatial  frequency  is  always 
attenuated  in  the  process  of  distortion,  a  safe  estimate  is  that 
the  spatial  frequencies  in  the  undistorted  image  must  have 
amplitudes  at  least  as  large  as  the  highest  amplitude  found  in  any 
of  the  individual  distortions.  It  may  also  be  observed  from 
studying  the  spectra  of  the  three  distortions  that  one  distortion 
may  produce  a  relatively  unattenuated  amplitude  for  one  spatial 
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frequency  whereas  another  spatial  frequency  in  that  same  distortion 
may  be  quite  heavily  attenuated.  It  might  therefore  be  expected 
that  by  selection  of  the  highest  amplitudes  for  all  spatial 
frequencies  from  all  distortions,  information  could  be  obtained 
which  would  be  superior  to  the  information  which  can  be  obtained 
from  any  one  distortion.  This  process  of  sorting  through  a 
number  of  distortions  and  selecting  the  highest  amplitudes  for 
each  spatial  frequency  has  been  termed  the  peak  spectrum 
selection  technique. 

Thirty  distortions  of  the  type  shown  in  Figs.  8  and  9  were 
constructed  and  their  spectra  generated.  The  spectra  were  then 
processed  by  the  peak  spectrum  selection  technique.  Figure  13 
shows  the  resultant  spectrum  obtained  by  selecting  the  peaks  or 
maximum  amplitudes  from  each  of  the  thirty  distortions.  It  may 
be  noted  that  the  spectrum  thus  obtained  is  a  very  close  match 
to  the  undistorted  spectrum,  the  envelope  of  which  is  shown  as  a 
dashed  line.  Phase  information  mast  also  be  processed.  Since  a 
symmetric  type  distortion  process  was  assumed,  that  is,  the  image 
is  equally  likely  shifted  left  or  right,  the  average  phase  shift 
is  0,  that  is,  it  is  equally  probable  that  the  phase  shift  will  be 
positive  or  negative.  A  simple  process  of  phase  averaging,  therefore, 
can  be  used  to  estimate  the  correct  phase  to  be  associated  with 
each  of  the  peak  amplitudes.  The  recovered  amplitude  and  phase 
information  may  now  be  utilized  to  construct  a  restoration  of  the 
distorted  image.  Figure  14  compares  the  image  restoration 
generated  in  this  matter  with  the  original  undistorted  image  shown 
with  the  dashed  lines.  While  the  recovery  obtained  was  not 
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complete,  the  reoorered  image  does  bear  oloee  resemblance  to  the 
undistorted  image  and  is  far  superior  to  any  one  of  the  thirty 
individual  distortions  utilised  in  the  process*  This  latter 
point  is  significant  since  it  indicates  that  this  prooess 
involved  something  more  than  simply  looking  through  a  large 
series  of  distortions  and  selecting  the  best  of  these  distorted 
images*  This  would  be  expected  from  the  fact  that  the  informa¬ 
tion  which  constitutes  the  recovered  image  was  obtained  not 
from  a  single  distortion  but  from  many  of  the  thirty  distortions* 

5.0  COMPUTER  PROGRAMING 

The  example  of  Section  4*0  clearly  indicated  the  value  of 
two-dimensional  Fourier  analysis  as  a  tool  in  the  processing 
of  distorted  images.  These  calculations  were  carried  out  on  a 
desk  calculator  and  required  a  great  deal  of  time  and  effort  in 
obtaining  the  recovered  image  of  Fig*  14.  It  was  recognized  that 
if  this  technique  was  to  have  application  to  two-dimensional 
images  having  a  large  number  of  picture  elements  and  for  a  large 
number  of  distortions,  the  hand  calculator  techniques  must  be 
replaced  by  computer  operations. 

The  Computer  Center  of  the  School  of  Science  and  Engineering, 
University  of  California,  San  Diego,  includes  a  Control  Data 
Corporation  1604  computer  with  a  160-A  satellite  computer* 
Peripheral  equipment  includes  an  IBM-088  card  reader ,  an  Analex 
high  speed  printer,  and  necessary  magnetic  tape  units.  Computer 
programs  for  this  project  have  been  written  in  Fortran  60*  The 
basio  computer  program  consists  of  taking  a  series  of  input  data 
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points  which  describe  the  image  at  points  on  an  XI  grid  and  per¬ 
forming  two-dimensional  xuurier  analysis  to  obtain  the  amplitude 
and  phase  of  each  spatial  frequency  comprising  the  image* 
Additional  subroutines  allow  for  operation  on  these  spatial 
frequencies  in  any  desired  fashion.  For  example  one  such  sub¬ 
routine  allows  the  application  of  the  peak  selection  technique 
and  the  phase  averaging  as  utilized  in  Section  4.0*  The  oomputer 
output  consists  of  tabulations  of  the  amplitude  and  phase  of 
each  of  the  spatial  frequencies  for  each  image  which  has  been 
inserted,  as  well  as  the  peak  amplitude  and  average  phases 
associated  with  restorations  achieved  by  the  peak  selection 
technique,  and  the  restored  image  obtained  by  taking  the  inverse 
transformation  of  this  peak  selection  phase  averaged  spectrum. 
Several  auxiliary  subroutines  for  displaying  the  output  ani 
input  information  have  been  developed  and  can  be  used  as  desired. 
As  an  example,  the  computer  can  take  the  output  image,  round 
each  picture  element  to  the  nearest  digit  between  one  and  ten, 
and  print  out  an  image  having  the  proper  XI  relationship  with 
this  ten  step  gray  scale. 

An  "image11  printout  was  developed  as  an  aid  to  immediate 
visual  evaluation  of  the  results.  Printed  patches  of  each 
character  available  in  the  Analex  printer  were  photometered. 

The  values  obtained  reflect  the  difference  in  inked  areas 
associated  with  the  different  characters.  The  maximum  reflect¬ 
ance  is,  of  course,  a  blank  and  the  minimum  reflectance  was 
found  to  be  the  character  £*  Characters  of  intermediate 
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reflectance  were  chosen  to  fill  out  a  ten  step  gray  ioale  as  follows* 


Character 

folfrvJ-Tf  P+awfonw 
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• 

0.95 
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0.84 
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0.81 

H 

0.78 

* 

0.75 

E 

0.72 

While  the  resulting  printed  images  are  lacking  in  contrast,  they 
do  allow  a  measure  of  visual  evaluation.  An  example  is  shown 
in  Fig.  15. 


6.0  MATHEMATICAL  MODEL  OF  ASTRONOMICAL  DISTORTIONS 


The  example  of  Section  4.0,  A  Geometric  Distortion  Model,  was 
extremely  valuable  as  a  means  of  investigating  the  application  of 
Fourier  techniques  to  restoration  processes.  It  does  not,  however, 
represent  an  accurate  model  of  the  distortion  process  which  would 
be  expected  by  a  large  aperture  optical  system.  The  geometric 
distortion  model  assumed  an  intensity  addition  of  images  formed 
by  the  various  sectors  of  the  entrance  pupil.  This  model,  therefore, 
neglects  the  effect  of  the  coherence  of  the  wave  front  over  the 
entire  entranoe  pupil.  While  it  is  perfectly  plausible  to  consider 
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Figure  15  Direct  Analex  printer  image  readout. 
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separately  the  images  formed  by  various  sectors  of  the  entrance 
pupil,  these  images  will  not  be  additive  in  intensity  beoause 
of  the  ooherenoe  of  the  electromagnetic  radiation*  A  more 
realistic  model  must  consider  the  addition  of  these  images  with 
due  consideration  to  the  relative  phase  of  the  electromagnetic 
radiation  as  received  from  the  different  sectors*  It  will  be 
shown  that  this  distinction  is  an  important  one. 

In  Appendix  II,  it  is  shown  that  any  two  points  on  the 
entrance  pupil  serve  to  generate  a  single  spatial  frequency 
which  has  a  sinusoidal  variation  in  a  direction  parallel  to  the 
vector  connecting  the  two  points.  For  example,  if  the  two 
points  on  the  entrance  pupil  were  chosen  to  fall  on  the  X  axis, 
the  resultant  image  would  be  described  by  the  equation 


F(x)  =  A-jAg  +  A-jAg 


(5) 


as  illustrated  in  Fig.  16.  The  terms  A^  and  A^  are  the 
amplitudes  of  the  electromagnetic  wave  front  at  the  two 
points  on  the  entrance  pupil,  D  is  the  distance  between  the 
two  points  on  the  entrance  pupil,  is  the  wave  length  of 
the  radiation,  f  is  the  focal  length  of  the  optical  system, 
x  is  the  spatial  dimension  in  the  image  plane,  and  A0  is  the 
phase  difference  between  the  electromagnetic  radiation  at 
the  two  points  on  the  entrance  pupil*  For  the  case  of 
atmospheric  turbulence  close  to  the  entrance  pupil, 

A^  and  A^  would  have  approximately  the  same  amplitude  as  in 
the  absence  of  turbulence,  and  the  primary  image  distortion 
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Figure  16  Intensity  distribution  from  two 
points  on  the  entrance  pupil. 
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vould  be  due  to  the  Nj&  term*  Appendix  II  goes  on  to  show  that 
the  total  contribution  to  a  given  spatial  frequency  may  be 
found  by  summing  all  of  the  sine  terms  of  the  type  shown 
in  Eq.  (5)  from  every  pair  of  points  on  the  entrance  pupil, 
separated  vectorially  by  the  distanoe  D(see  Eqs.  11-15  and 
11—16)  •  Figure  17  illustrates  this  point  by  showing  that  area 
of  the  entrance  pupil  which  is  responsible  for  the  generation 
of  a  spatial  frequency  D  *  This  has  been  termed  a  spatial 

X 

frequency  zone. 

As  in  the  case  of  the  geometric  distortion,  a  maximum 
amplitude  is  obtained  when  all  of  the  A0*s  are  equal  or  equal 
to  0.  This  is  the  case  of  an  undistorted  image  where  the 
electromagnetic  wave  front  is  still  plane  at  the  time  it 
reaches  the  entrance  pupil.  Thus  it  is  once  again  apparent 
that  the  effect  of  atmospheric  turbulence  will  be  to  attenuate 
the  amplitude  of  any  given  spatial  frequency.  The  peak 
spectrum  selection  technique  utilized  in  restoration  of  the 
image  distortions  in  Section  4.0  is  therefore  applicable  to 
this  more  refined  model  of  distortions  due  to  atmospheric 
turbulence. 

The  distortion  process  may  be  visualized  with  reference 
to  Fig*  16  by  imagining  a  phase  distortion  map  superimposed 
over  the  entrance  pupil  of  the  optical  system.  The  amplitude 
of  any  spatial  frequency  in  the  resulting  distorted  image  is 
found  by  summing  all  of  the  sine  waves  of  Eq.  (5)  for  every 
pair  of  points  separated  by  a  distance  D  as  defined  by  the 
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zones  of  Fig.  17.  A  power  spectrum  m ay  be  associated  with  the  phase 
distortion  map.^  This  power  spectrum  will  fall  off  as  frequency 
is  increased.  The  significance  of  this  fact  is  that  there  will  be 
no  discontinuities  or  very  rapid  changes  in  the  phase  map  over 
short  distances.  Two  points  on  the  phase  map  separated  by  a  short 
distance  will,  therefore,  be  highly  correlated.  Since  the  very 
low  frequencies  passed  by  the  optical  system  correspond  to  pairs 
of  points  separated  by  very  short  distance,  there  will  be  a 
tendency  for  the  very  low  frequencies  to  be  passed  with  small 
attenuation.  For  the  higher  frequencies  which  relate  to  points 
separated  by  increased  distance,  there  will  be  an  increased 
probability  of  having  phase  difference  between  the  two  points  and, 
therefore,  it  is  to  be  expected  that  as  the  frequency  increases, 
the  attenuation  will  increase.  The  increase  in  attenuation  with 
frequency  will  continue  until  a  frequency  is  reached  at  which  the 
separation  between  the  pairs  of  points  which  form  this  spatial 
frequency  is  sufficiently  large,  such  that  the  phase  variation 
at  the  two  points  becomes  independent.  From  this  point  on 
there  will  be  no  increase  in  attenuation  with  frequency. 

There  is  a  very  great  significance  to  the  assertion  of  the 
previous  paragraph.  Previous  geometric  models  of  the  distortion 
process  have  suggested  that  the  image  of  a  point  source  in  the 
presence  of  turbulence  would  be  normally  distributed  in  two 
dimensions.  The  Fourier  transformation  of  a  two-dimensional 
normal  distribution  is  also  a  two-dimensional  normal  distribution. 


•^S.  H.  Reiger,  "Seeing  through  the  atmosphere,*  Proceedings  of  a 

Symposium,  December  1962,  Rand  Corporation  Memorandum  RM-3294-PR* 
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Such  a  model  suggests  that  there  is  an  increase  of  attenuation  with 
frequency  without  limit*  This  type  of  model  would,  therefore, 
impose  a  practical  limit  to  the  highest  frequencies  which  could 
ever  hope  to  be  recovered  or  restored.  The  fact  that  the 
more  realistic  model  indicates  a  leveling  off  of  attenuation 
with  the  higher  frequencies  suggests  that  no  such  limit  does, 
in  fact,  exist. 

In  addition  to  offering  an  explanation  of  the  distortion 
process  in  the  presence  of  atmospheric  turbulence,  the 
derivations  of  Appendix  II  also  offer  a  ready  means  of 
determining  the  transfer  function  of  a  diffraction  limited 
optical  system.  Lince  the  contribution  to  a  given  spatial 
frequency  can  be  associated  with  zones  on  the  entrance  pupil 
bounding  all  points  separated  by  distance  D  =  f A, 
as  illustrated  in  Fig.  17,  the  defining  of  these  zones  and 
the  measurement  of  this  area  relative  to  the  totil  area  of 
the  entrance  pupil  is,  directly,  the  transmittance  of  the 
optical  system  for  this  given  spatial  frequency.  This  concept 
of  zones  associated  with  non-coherent  optical  systems  is 
suggestive  that  where  primary  concern  is  for  objects  of  high 
spatial  frequency,  techniques  of  apodization,  that  is, 
masking  of  the  entrance  pupil,  may  be  effective  in  increasing 
the  amplitude  of  the  spatial  frequencies  relative  to  the 
average  ambient  amplitude. 


7.0  REAL  TURBULENCE  DISTORTIONS  AND  RESTORATIONS 


The  peak  spectrum  selection  technique  introduced  in  Section  4.0 
has  also  been  utilized  in  attempts  to  restore  images  distorted 
by  real  turbulence.  A  six-inch  diameter  astronomical  objective  was 
mounted  so  as  to  form  an  image  in  the  plane  of  a  mechanically 
driven  aperture.  The  flux  passing  through  the  aperture  was 
collected  on  the  photooathode  of  a  multiplier  phototube  whose 
output  was  displayed  on  an  oscilloscope.  The  scanning  aperture 
performed  a  single  horizontal  scan  and  the  resulting  oscilloscope 
waveform  was  photographically  recorded.  The  scanning  apparatus 
was  equipment  available  at  the  Laboratory  as  a  result  of  work 
on  a  prior  photoelectric  scanning  problem.  A  hot  plate  was 
mounted  in  front  of  and  below  the  objective  lens  to  produce 
turbulence  close  to  the  entrance  pupil.  The  single  horizontal 
‘  scan  was  sufficiently  fast  so  that  the  image  was  essentially 
fixed  during  the  entire  scan.  In  order  to  obtain  high  flux 
levels  and  therefore  minimize  the  noise  level  relative  to  the 
signal  level,  the  object  to  be  imaged  was  chosen  to  be  a 
filament  lamp.  Since  the  scan  was  performed  in  one  dimension 
only,  it  was  necessary  that  the  image  have  flux  variation  in 
only  one  dimension.  The  image  of  the  filament  lamp  satisfied 
this  requirement  in  that  for  a  considerable  region  in  the 
image  plane,  it  had  the  general  appearance  of  being  three 
vertical  bars.  Figure  18  shows  thirty  osoillosoope 
photographs,  each  representing  one  horizontal  trace  through  the 
image  of  the  filament  lamp.  The  fourth  recording  in  this 
series  was  taken  with  the  hot  plate  removed  in  order  to  show 
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Figures  18  through  22 
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the  nature  of  the  undistorted  inage*  The  reduction  of  information 
content  by  the  atmospheric  turbulence  can  be  visualised  by 
imagining  that  the  essential  information  which  was  desired  from 
the  image  was  how  many  peaks  are  present  in  the  image,  and 
secondly  what  is  the  relative  order  of  amplitude  of  the  three 
peaks,  i.e.,  which  is  largest,  which  is  next  largest,  and  whieh 
is  smallest* 

Before  attempting  the  peak  spectrum  selection  technique, 
some  more  common  forms  of  processing  were  performed.  For 
example,  Fig.  19  shows  the  addition  or  integration  of  the 
twenty-nine  distorted  waveforms.  The  result  is  an  image  in 
which  the  third  peak  is  completely  missing*  Thus,  in  terms 
of  the  imagined  information  desired  from  the  image,  this  type 
of  processing  is  unsuccessful. 

One  component  of  a  turbulence  distortion  can  be  visualized 
as  a  gross  shift  of  the  image.  It  would,  therefore,  seem 
reasonable  that  this  shift  could  cause  a  considerable  loss  of 
information  in  the  process  of  direct  integration  previously 
described*  For  this  reason  a  second  integration  was  performed, 
but  this  time  the  center  of  gravity,  i*e.,  the  first  moment  of 
the  flux  distribution  was  determined  and  all  images  were  shifted 
to  have  a  common  center  of  gravity,  prior  to  integration* 

Figure  20  shows  a  result  of  this  second  integration*  There  is 
no  significant  improvement  associated  with  this  form  of  process¬ 
ing. 

The  twenty-nine  distortions  were  then  subjected  to  Fourier 
analysis  with  the  aid  of  the  GDC  1604  computer  and  the  peak 
spectrum  selection  technique  was  applied.  This  was  first 
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accomplished  without  using  the  center  of  gravity  shift  and  the 
results  are  shown  in  Fig,  21,  The  third  peak  is  now  present}  the 
three  peaks  have  approximately  correct  relative  amplitudes 
although  considerable  distortion  still  exists.  Figure  22  shows 
the  results  of  applying  the  peak  spectrum  selection  technique 
where  the  phase  estimates  were  made  after  all  images  were 
shifted  to  common  center  of  gravity.  Comparison  with  Fig,  18 
shows  that  while  the  restoration  of  detail  is  not  complete,  i.e., 
it  is  lacking  in  contrast,  much  of  the  detail  of  the  original 
image  has  been  recovered. 


8.0  TIME  INVARIANT  DISTORTIONS 

8.1  Techniques  and  Noise  Limitations 

The  preceding  section  dealt  with  the  case  of  time  variant 
image  distortion,  i.e,,  the  situation  in  which  a  time  sequence 
of  recorded  images  shows  a  different  distortion  in  each 
exposure.  Another  case  of  interest  is  the  time  invariant  dis¬ 
tortion  in  which  a  time  sequence  of  image  recordings  shows  an 
identical  distortion  for  each  exposure.  Lens  aberrations  would 
be  of  this  class  of  distortion.  Atmospheric  turbulence  oan 
also  produce  this  type  of  distortion,  for  example,  in  the  case 
where  a  long  time  exposure  is  required.  In  a  long  time 
exposure  a  large  statistical  sample  of  the  turbulence  is 
averaged  and  there  is  therefore  little  variability  from  image 
recording  to  image  recording.  This  type  of  distortion  is 
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therefore  characterized  by  a  transfer  function  which  is  temporally 
invariant.  Figure  23  shows  a  hypothetical  transfer  function 
for  this  type  of  distortion. 

If  the  transfer  function  is  known,  restoration  can  be 
accomplished  by  amplifying  each  of  the  spatial  frequencies  by  an 
amount  which  exactly  compensates  for  the  attenuation  which  it 
has  suffered.  The  limitation  in  quality  for  this  type  of 
restoration  lies  in  the  presence  of  noise.  Any  real  sensor 
recording  an  image  is  ultimately  limited  by  some  type  of  noise. 

In  the  example  shown  in  Fig.  23,  the  noise  is  assumed  to  be 
white  or  flat  with  frequency.  In  the  process  of  amplifying  each 
of  the  spatial  frequencies  to  compensate  for  its  attenuation, 
the  amplitude  of  the  noise  of  that  spatial  frequency  will  also 
be  amrvLified  by  an  identical  amount.  Figure  23  shows  the 

•is 

noise  spectrum  which  would  result  after  compensation  for  the 
transfer  function  of  that  example.  A  spatial  frequency  which 
had  an  amplitude  equal  to  the  noise  prior  to  restoration  will 
have  an  amplitude  equal  to  noise  after  restoration.  It  is 
conceptually  important  that  the  information  content  of  a  single 
image  is  not  increased  by  any  such  processing.  However,  where 
simple  additive  noise  is  the  limitation  in  the  restoration, 
the  signal-to-noise  ratio  can  be  increased  by  integration 
processes  well  established  in  the  fields  of  radar  and  communica¬ 


tion, 
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Figure  23  Typical  transfer  function. 
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8.2  Processing  of  Noisy  Inages 

A  demonstration  of  the  improvement  of  sign&l-to-noiee  ratio  by 
integration  was  performed.  An  image  orthicon  television  chain  was 
used  to  view  a  test  pattern.  The  chain  was  purposely  adjusted 
to  give  an  extremely  poor  signal-to-noise  ratio  on  the  video 
monitor.  The  upper  left  photograph  in  Fig.  24  shows  the 
photographic  recording  of  a  monitor  image  of  a  single  frame  (l/30 
second)  of  the  test  pattern.  As  can  be  seen,  the  noise  level  was 
sufficiently  high  so  that  little,  if  any,  information  is  present 
in  this  recorded  image.  Integration  was  then  accomplished  by 
increasing  the  exposure  time  of  the  camera  used  to  photograph  the 
monitor.  For  example,  the  second  image  shows  an  exposure  time 
four  times  that  of  the  first  image  with  the  fi  number  of  the 
optical  system  adjusted  to  maintain  constant  integrated  flux.  The 
sequence  of  pictures  demonstrates  the  effect  of  increasing  the 
number  of  frames  which  are  averaged  or  integrated.  Figure  24 
shows  clearly  the  increase  in  information  content  which  can  be 
obtained  by  increasing  the  period  of  observation  or  the  number 
of  images  which  are  processed.  If  it  is  imagined  that  the  upper 
left-hand  picture  representing  an  exposure  of  one  frame  was  the 
image  which  resulted  from  the  correction  of  a  temporally 
invariant  distortion,  then  the  pictures  which  follow  indicate  the 
improvement  in  image  information  which  could  be  obtained  as  a 
result  of  processing  of  multiple  images. 


64 
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Figure  24  Illustration  of  the  effect  of  time  averaging.  The 
original  noise  level  is  128  times  normal  image 
orthicon  noise  and  the  number  below  each 
picture  indicates  the  number  of  frames  which 
have  been  integrated. 
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8.3  Example  of  a  Tine  Invariant  Restoration 

The  techniques  involved  in  the  restoration  of  temporally 
invariant  distortions  were  studied  by  means  of  simple  mathematical 
distortions  of  this  class.  One  example  of  a  temporally  invariant 
distortion  is  the  case  of  a  photograph  taken  while  the  object  to 
be  photographed  is  in  motion  all  during  the  time  of  the  exposure. 
This  would  be  the  equivalent,  for  example,  of  a  reconnaissance 
photograph  taken  under  conditions  of  inadequate  image  motion 
compensation.  Figiire  25  shows  a  distortion  of  this  type  which 
was  mathematically  constructed.  Here  the  block  letters  ARPA 
were  smeared  horizontally  by  an  amount  which  corresponds  to  an 
image  motion  during  the  exposure  equal  to  1  l/2  times  the  width 
of  a  block  letter.  The  resulting  image  distortion  is  shown  in  the 
top  of  Fig.  25.  This  distorted  image  along  with  information  as 
to  the  nature  of  the  distortion  process  were  fed  into  the 
computer  which  then  performed  Fourier  analysis  of  the  distorted 
image,  amplification  of  each  of  the  Fourier  components  to 
compensate  for  their  attenuation,  and  the  construction  of  the 
restored  image  which  is  shown  in  the  bottom  of  Fig.  25.  In  this 
"noiseless  case"  the  image  recovery  was  accurate  to  the  eighth 
significant  figure  and  was  limited  only  by  the  extent  to  which 
the  finite  Fourier  series  could  approximate  the  image  in 
question. 
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Figure  25  image  distortion  and  recovery. 
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9.0  THE  FILM  SCANNER 

The  preceding  sections  hare  described  "techniques  of  image 
restoration  and  shown  examples  of  distortions  and  restorations 
by  means  of  which  these  techniques  have  been  explored  and  proved. 
It  was  recognized  at  an  early  state  in  the  project  that  the 
restoration  of  atmospheric-turbulence-distorted  images  by  these 
techniques  would  require  some  type  of  automatic  accumulation  of 
input  data  for  the  computer.  The  most  easily  adapted  sensor  for 
recording  the  distorted  images  is  photographic  film.  For  this 
reason  it  was  decided  that  the  initial  trials  of  real  image 
restorations  would  be  accomplished  from  images  recorded  on 
photographic  film.  The  number  of  resolution  elements  associated 
with  the  image  distortions  virtually  prohibits  the  hand 
processing  of  these  films  in  order  to  obtain  input  data  for  the 
computer.  A  considerable  projeot  effort  has  therefore  been 
directed  toward  the  construction  of  a  film  scanner  capable  of 
automatically  reducing  the  input  data  into  computer  format. 

Figure  26  is  a  photograph  of  the  completed  film  scanner  and 
associated  electronics.  At  the  left  rear  of  the  picture  may  be 
seen  a  projection  system  through  which  16  mm  film  is  fed.  Films 
of  35  mm  or  70  mm  would  be  printed  on  16  mm  format  for  use  in 
this  scanner.  Scanner  heads  to  accommodate  the  larger  formats 
oan  be  constructed  if  required.  The  image  is  projected  onto 
the  rectangular,  flat,  white  metal  mask  (or  screen)  shown  in  the 
middle  of  the  picture.  The  mask  contains  a  small  aperture  and 
the  flux  passing  through  the  aperture  is  Imaged  onto  the 


Figure  26  Film  scanner  (in  rear)  and  associated  electronics. 
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photooathode  of  a  multiplier  phototube  which  is  mounted  in  the  black 
box  behind  the  white  mask. 

Figure  27  shows  a  close-up  of  the  projection  system.  The 
operation  of  the  system  will  be  explained  with  the  aid  of  the 
block  diagram  shown  in  Fig.  28.  The  scanning  is  accomplished  by 
physical  movement  of  the  film  in  the  image  plane.  A  horizontal 
scan  is  accomplished  in  a  series  of  discrete  steps.  Stepping 
pulses  are  supplied  by  the  motor  driver  unit  and  activate 
stepping  motors  in  the  projection  unit.  The  horizontal  counter 
counts  the  number  of  horizontal  steps  which  have  been  taken. 

When  the  desired  number  of  steps  have  been  accomplished  the 
motor  driver  unit  supplies  a  series  of  rapid  pulses  which  cause 
the  stepping  motor  to  retrace  to  its  original  position,  and  at  the 
same  time  take  one  vertical  step.  The  horizontal  scan  is  then 
repeated.  When  the  desired  number  of  vertical  steps  have  been 
completed,  the  motor  driver  supplies  a  series  of  rapid  pulses 
to  the  vertical  motor  which  cause  it  to  retrace  to  its  starting 
position.  The  system  is  now  ready  for  a  new  film  frame. 

The  output  of  the  nrultiplier  phototube  is  fed  into  a  voltage- 
to-frequency  converter.  The  converter  output  is  then  counted.  The 
digital  output  of  the  counter  is  serialized  so  that  it  may  be  fed 
directly  to  an  IEM  card  punch.  The  system  in  this  way  accomplishes 
automatic  recording  on  I  EM  cards  of  the  images  contained  on  the 
film. 

Where  commercial  units  were  employed,  their  identifying 
trade  names  and  serial  numbers  are  shown  on  the  block  diagram.  The 
electronic  circuits  designed  and  constructed  by  the  Visibility  Labora¬ 
tory  are  shown  in  Figs.  29  through  36. 


Figure  27  Close-up  of  film  scanner  head. 
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Figure  28  Block  diagram  scanner  system. 
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Figure  29  Circuit  Diagram  -  Serializer. 


Figure  30  Circuit  Diagram  -  BCD  Converter. 
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Figure  31  Circuit  Diagram  —  Reset  Control. 


Figure  32  Circuit  Diagram  -  Forward  8  Reverse  Pulse  Generator. 
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Figure  33  Circuit  Diagram  -  Horizontal  &  Vertical  Motor  Control. 
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Figure  35  Circuit  Diagram  -  Switch. 
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Figure  36  Circuit  Diagram  -  Horizontal  8  Vertical  Counter*. 
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At  the  time  of  writing  of  this  report,  the  oonetruotion  and 
assembly  of  the  film  scanner  system  has  been  completed  and  effort 
is  presently  directed  toward  the  scanning  of  test  films  used  to 
check  the  reliability  of  the  system. 

10.0  FUTURE  PLANS 

The  computer  program  used  to  accomplish  the  restoration  of 
the  smeared  ARPA  image  has  been  rewritten  to  allow  for  smear  in 
two  dimensions.  A  mathematical  distortion  of  the  block  letters, 
ARPA,  corresponding  to  a  two-dimensional  smear  will  be  fed  into 
the  computer  as  a  test  of  this  program.  At  the  same  time, 
photographically  recorded  two-dimensional  smears  of  the  block 
letters  ARPA  are  being  prepared.  These  are  accomplished  by 
masking  the  entrance  pupil  of  the  recording  camera  to  the  form  of 
a  square  and  then  defocusing  the  image  so  that  every  point  in 
object  space  becomes  a  square  in  the  image  plane.  As  in  the 
case  of  the  one-dimensional  smear,  the  smear  function  will  be 
sufficiently  large,  such  that  the  block  letters  ARPA  cannot  be 
identified  in  the  distorted  images. 

Because  of  the  presence  of  film  noise,  the  restoration  of  & 
single  distortion  from  a  single  film  frame  may  result  in  an 
image  of  very  high  noise  level.  For  this  reason  multiple  images 
will  be  recorded,  that  is,  a  sequence  of  film  frames  are  being 
prepared,  each  of  which  is  an  identical  smear  of  the  block  letters 
ARPA.  The  computer  program  will  therefore  involve  not  only 
restoration  of  each  frame  but  also  an  integration  of  the  resulting 
recovered  images  to  form  a  single  image  of  high  eignal-to-neiee 
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When  the  temporally  invariant  distortions  of  this  type  have 
been  successfully  recovered!  the  next  step  will  be  to  generate 
film  strips  with  the  optical  system  used  to  produce  atmospherically 
distorted  images  in  the  laboratory.  The  flux  levels  and  time 
exposures  will  be  chosen  so  as  to  give  large  statistical  sampling 
of  the  turbulence  and  therefore  produce  temporally  invariant 
distortions.  After  these  distortions  have  been  successfully 
processed,  additional  atmospheric  distortions  will  be  generated 
of  the  temporally  variant  variety.  These  film  strips  will  be 
scanned  and  fed  into  the  computer  and  the  two-dimensional 
version  of  the  peak  spectrum  selection  technique  will  be  used  to 
attempt  restorations. 

The  above  trials  will  represent  a  checkout  of  the  entire 
system  from  the  film  scanner  through  the  computer  programing  and 
will  constitute  a  test  of  the  restoration  techniques.  When  these 
techniques  and  the  system  have  been  demonstrated  to  be  performing 
adequately,  the  system  will  be  ready  for  attempted  restorations 
of  some  real  space  photographs.  While  the  system  as  presently 
constructed  is  inadequate  for  any  volume  processing,  it  should  be 
extremely  adequate  for  demonstrating  that  linage  improvement  can 
be  obtained  by  the  techniques  which  have  been  described.  The  task 
of  engineering  a  system  capable  of  operation  with  an  astronomical 
system  should  not  be  initiated  until  the  simple  system  which  has 
been  constructed  has  clearly  demonstrated  the  feasibility  of 
accomplishing  this  type  of  image  restoration* 
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'fov'vvMfltwi  $  fattwi  sgv’^s  oy k  a  two  -  dWv^mitvml  sfT  of 
dUorrte  ptfvrtT  valves. 

As  %y»  Sttfiw  (s  Ms  cUViVfcTiw  \%  lo*j9i  o*  HilcUbfftY^ . 

Ito  dtfwwMtfiiY*  u )‘»U  b?  1<(T*y  »*  tvit 

cU  rwocT\ov\ : 

Kv>o\m\a: 

AH 

2-iiO  Ci>S  1TCW  Cos  Z  ire  YV  =  0  J  W\  1 V)  (l.|) 

a  a 

-  A;Y«  =  n^o,iL  (1*2) 

2.  ^ 

-  ajY n=n  =  o,&  64) 
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a-i 

Ic»#  Sw  Siv»2tun  «  Oj  vn*w  fa.i) 

0>  A 

*A;W*n/Oor^.  (14) 

*  Oj  wi«\n-*OoviL  (i.i) 

1 

a-i 

2t.0  s»v)ixi£Q  tosinin»0  (5) 

<k  & 


4 1*  -  _|?  Lx.  .  bO,P 
P 


©te  *  ^  TT^fc  -  ZTrit 

L*  P 


1>r»y»«:  1  L>t_  .  1  *  0> CJ 


4>1  *  «  ZTT  1 


W**t«  Lx  is “tW  Unoty  o(  -fV»c  p*<iod  of  1Vi<  function  on  *  And 
L^  is  IVk  Un^V)  of  IV  period  on  u. 
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p  a*  A  <xv«  iV»e  vtvwttaY*  oV  i\AA»^v»AonT  po\vtr^  on  \  awA  ^ 
oCT  wWvtU  1W  CiMftdflw  is  dtfivwA .  Total  points  ow  %  anA  c|  <w 
(p+i)  av»A 

EKpfcnA'm^  $  ify)  U)ttU  y«j1?(TTo  0: 

? (6,40  *  ^  it0  C&l(ti  cos ie  +  'SifdOsWtie] 


£xp*nA\v\^  A  c  (&)  &nA  (4>V* 

"  Zjlo  [Cij  COS j<t>  -1-  T>ij  siv>  j  <t>] 

*  ^y-o  [EildSit  *  TljSiv>]<p] 
?(8v,^iV-Zt'6  Zj^  (Cij  cosL©k  cos]<5>!  t 
■‘■‘Pi]  Cost  Oil  Sm^L  +  £(,;  Sini  6h  CDS  j  <£i  t 
+  Tij  svniQk  Sinj^!^  (H) 

To  sViiw  dWl  dVis  fiAnctioirt  YfyY«s«mrfs  tVi*  vmimkIu  orient *d 
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s\v»  and  cos  <oav«fa\w* ,  IW  function  can  b<  wf* rttm  .* 

f  (©k>iV-  \  2. if.  Zjt  [^li-^^Oek+j'fi)1’ 

*  (x>i]t  E;j)s\v«(iek  -+  j4>iV 

+•  (-*Di j  +  E  ^)s'ivi  ( c  e  it  -  ]  4>i)  J 


To  dtfwvmvu  C  ij  VAUH'ipl^  ( *f)  b>  ^  (cos  rn  ©fc  COS  X)  4>  i') 

Avid  SUVA  OUW  1W  p  Mid  Cj  ivt<Up«V»A*nt  OC  OlVIfl  C|  dftjs 
powits: 

p*l  q  - 1  .  * 

Zk.0  2 1*,  f(©k,^i)  cos  w6k  Cosn<))i  * 

*  Z  it .  2  jt  2  kM  £  1-0  Cy  COS  (0k  COS  Kl0kCOS jK  <OS«l*j. 


h\\  oiWtlmws  vanish  sivut  Hmj  Contain  a  summation  of 


sin- Cos  f»ms  anA  -TV»t  sumtudfion  yttuMfc  in  0,  from  (*3} . 
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Tbt  about  «<|uortittvi  ba*  valu«  ovil^  u)bm  i.  -  vr»  av>d  j«w 
swt»ulTavwovisli| ,  frow  ( I  .f) 

2k£i  cow  Oh  Coijta  = 

=  Cos'c6h 

Ci^  S  I.lt^l  P  COS  (9  k  cosi^x  ^  Pc 

Ik  Zi  cosMek  tos'j^x  Gc 


WW  i  ?  09_|.  av>A  j  ^  0,^.  ,  fot>«\(\,i)i 

0  c  -  -£4-  av»<\  C*:  ;  =  iL  y, 

h  ^  PS 

bib***  i  *  0  ot^.  Q\rt<l  j  ^  0  ov  ^ 

OY  i  ^0  0¥  avt/l  j  =  0  OY  4  ,  (Vottl 


=  -P-3-  ^  i\  ~  _2i_  P% 

z  4  PT 

W*IW  1  =  0  OY^  <XV»A  j  s  0  6Y  |  ,  pYOVrt  6*i) 


Gc  =  PS  a*d  Cis  -  J_  Fc 

ri  4  PS 


To  <w«\ftpu(4)  bu  cos  W0fe  si«  ovti 


Suw. 


1.  f  (§*,4*0  Coswi^h  $wi  n^i  - 
"  2-k  (josL©^  CoswOfcsWijcJ^sivin^i 

U)i1V>  vaU«  ioViqy*  i:hi  <*v\A  j  -  n  ^  0>3_  f  yowi  (2*2)5 

g  2 

T>ij  -  itii  £[e^t)  cost  BhSlttj^i  c  Fa 
Z  W  X 1  COS1  c  ©H  SivV1  j  ^  el 


U)Vk9Yl 


i  *  Ooy|.  av\a  ]  ^  0,j|.  ,fY0*>(|.z),(2.2) 


<3d  r  omA  'T)  t  :  =  JL  F<i 

4  J  pc, 

U)V\M  L  2  0  OY  omA  j  4  fl  OY  | 
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To<WTw«nfl<  Ei/j  ItawUipl^  (q)  bi^  sirtv*0b  dosn^L  am!  5UM: 

£k£i  P(6k,$i)siv» w»ew  cosn^i  * 

=  Zli^i  E.iiZx  E^-smieit  S»v>vn0it  C03j^l  C0*n<j>! 

UjVWa  valu<  u)W<  i.  -  no  «*  0  ov.|  and  \  =  n  ,  from  (2.2) 

£1;  8  Z  >tX  1  sime*  <0S}4>1  _  Yl 

Z  it  X 1  s\vP  i.  ek  cos1]^!  G< 

U)Vwv>  i  ■*  OOv^.  Qv>d  j  /  0  ov|_  j  from  (l.i) ,  (l.l) 

G<  =  pi  and  £:;  »  4.  F« 

"  ■»  pci 

WU»vi  i  *  0  Or  J>  and  ]  *  0  OVjJ.  ,  Pxo^('|.^)>(2.l) 

G«c_P_4  avid  -  2_ 

J  p| 
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To  dtfwwwt  bi^  s »v%  w\  0  s’iy>  owd 


^SU tf\ : 


Zk  Z.J.  f (6^ (4>i)  s'k»  wi ©k  s’m  y\ * 

~Zl  £j  ^wZ.i^k,4>iV'^^w0ilSm^lSivjn<(> 
VxiVtv\  U)Vwy«  i  '  wt  *  0,_£  avwJi  j  -  ^  ^  0  , 

(Vftv*  (2 .2") ; 

Fj,j  =  Z w  Z.  x  f  (6k  /O  s>o  t  Q w  ,  Tf 

ZkZj.  SWfci©iiS\**]<t>i  <3* 

L0^er\  l  -  Yn  ■*  0  ov  £  Qwd  j  -  n  ^  0  6Y  i|_  ,(yow'  : 
Gf  ff  P4  o^d  T:;  -  Ff 

*  J  pcj  r 


Wuliow  WwlUu^ftiVwwUQO 
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Appendix  II. 

si  gflMUaaft  Isl  treagfor  Xtoailaa  si  &  fllffrwttQB 

limited  optical  system. 

Assume  that  a  diffraction  limited  optical  system  is  Imaging 
an  infinitely  distant*  monochromatic  point  source.  The  plane 
electromagnetic  wave  incident  on  the  entrance  pupil*  in  the 
absence  of  turbulence,  can  be  described  by  the  equation 

a.  =  &  s\vi  ( t  +•  4>'') . 


Suppose,  however*  that  in  the  neighborhood  of  the  entrance  pupil, 
the  plane  wave  passed  through  a  thin  lamina  having  a  spatially 
varying  index  of  refraction.  The  emerging  wavefront  would  then 
have  the  form 


where  x  and  y  are  the  horizontal  and  vertical  directions  in  the 
plane  of  the  entrance  pupil.  This  phase  disturbance  will 
result  in  an  image  distortion*  a  mathematical  derivation  of  which 
will  now  be  developed. 

The  nature  of  the  image  can  be  determined  by  the  application 
of  Huygen's  principle  in  which  each  point  on  the  wavefront  will 
be  considered  to  be  an  isotropic  radiator.  The  propagation  in 
a  direction  (oC  *^)  will  then  be  examined*  where  &  and  Q 
are  angles  measured  with  respect  to  the  x  and  y  axes  respectively. 


(II-l) 


(II-2) 
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The  optical  system  would  produce  an  intensity  map  as  a  function 
of  OC  and  such  that  the  intensity  is  proportional 

to  the  square  of  the  amplitude  of  the  vector  sum  over  the 
wavefront  in  a  plane  inclined  at  the  angle  (be,  .  The 
constant  of  proportionality  involves  the  geometry  of  the  optical 
system  and  is  not  relevant  to  this  particular  derivation. 

For  the  application  of  interest,  the  objects  to  be  imaged 
have  angular  subtense  on  the  order  of  seconds  of  arc,  and  the 
sine  and  tangent  functions  may  be  represented  by  the  angle 
expressed  in  radians.  Using  the  central  axis  as  a  reference, 
a  point  x,y  on  the  wavefront  will  have  a  path  length  increase 
in  the  inclined  plane  relative  to  the  point  0,0  of 


A  d  =  Y  (II— 3) 

This  increased  path  length  results  in  a  phase  shift  relative 
to  the  entrance  pupil  of 


y  -  oi'v  ±  (n-4) 

A. 

Equation  (II— 2)  can  therefore  be  rewritten  for  the  case  of 
propagation  in  the  direction (o(,  ^>)  as, 

siri[(Dt  (II“5) 

The  square  of  the  amplitude  of  the  vector  sum  of  Eq.  (II-5) 

over  all  x  and  y  within  the  limits  of  the  entrance  pupil  will 

be  proportional  to  the  intensity  of  the  image  at  the  point 
corresponding  to  the  direction 
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The  Image  plane  intensity  ia  therefore 


* 


+ 


(II-6J 


v. 


where  A(x,y)  ia  zero  in  the  region  of  x,y  outaide  the  boundary 
of  the  entrance  pupil.  By  recognizing  the  concept  of  dummy 
variables  of  integration,  (II-6)  can  be  rewritten  as 


I=  1.1  • 


■  c.os jaw  +s^.  +4>(w,i')  dwik  *■ 

+  1-.  +  +  <W^\  d*dp 


(n-7) 


o  oo  (? 00  r  l 
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Since  the  integration  with  respect  to  x  and  j  ie  Independent  of 
the  Integration  with  respect  to  w  and  z,  Eq»  (Hr?)  nay  be 
rewritten  as 


iZi*  m*.  +4’!^] 

•  Cosj^ +  $£  +  -v 


■ins  !•  S-*  *(***W[«UL  + 

si« 4* 


(H-8) 


By  using  the  trigonometric  identity 

cos  a  c&s  b  f  sin  a  s'mb  *•  Cos  ( a.  - 1>)  (n-9) 


Eq,  (II-8)  can  be  rewritten  as 

1  *  O-t  it  il  A  aMcos]**  tty  + 


-  2^£L  -^L  n d w c\ i 


(11-10) 


Let 

w*  ot  +  e^  ,  dw  * 

)  d  *  <•  d€ij 

i«l«L  I-®  5-®  ’^  + 


(ii-ii) 
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By  using  the  trigonometric  identity 

cosfa+b'}  »  COSO  co$b  -atoiAzinb 


(11-12) 


Eq.  (11-11)  becomes 
r  *  »“ 

r 


:0-»  cos[«es  (11-13) 

•  cosjj^Cocjiy)  -  +  €i|)]  <h<At|"j  d«»d€i|  + 

+  LLsiy,[^ 

•  -4>  (\+6x,t|+€^]  dtdljl  <U%  d*y 


The  terns 


Gr,(*,(j>)  *  cost^f*  *  H*"]  ] 


St(*,£)  -SmQ2^  +  ^jW  J 

may  be  recognized  ae  spatial  frequency  components  where 
and  the  double  integral  terms  which  multiply  them  in  Eq.  ( 11-13) 
is  the  transmittance  associated  with  the  spatial  frequency 
&.«iv  The  transmittance  of  the  cosine  term  is 

Tc  *  $-*  i-«  A(*  +  4*,ll  +  €i|V 

•  (oc -<K*+e*)^+  e(j)]d^du 


(11-14) 


(11-15) 
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(11-16) 


and  for  the  sine  term 

T*  -  §-*  • 

It  must  be  remembered  that  A  has  nonzero  value  only  over  the 
region  of  the  entrance  pupil.  The  transmittance  may  be  seen 
to  be  determined  by  pairs  of  points  on  the  entrance  pupil 
separated  by  a  distance  6  %  in  the  x  direction,  and  £  ij 
in  the  y  direction.  It  Is  further  noted  that  It  is  the 
product  of  the  amplitudes  of  the  pairs  of  points  and  the 
difference  in  electromagnetic  phase  between  the  pairs  of  points 
which  determine  the  transmittance. 

This  leads  directly  to  the  concept  of  Spatial  Freouencv 
Zones  on  the  entrance  pupil.  For  a  given  spatial  frequency 
(W},  those  regions  of  the  entrance  pupil  which  contribute  to 
the  generation  of  this  spatial  frequency  are  defined  by  bounding 
all  pairs  of  points  which  are  separated  by  a  distance  £%  in 
the  x  direction  and  in  the  y  direction. 

Equations  ( 11-15)  and  (11-16)  are  the  transmittance  values 
which  result  from  both  the  atmospheric  turbulence  and  the 
diffraction  limit  of  the  optical  system.  To  make  this  point  clear, 
the  transmittance  in  the  absence  of  the  phase  disturbance  will  be 
derived.  When  h$  =  0 

Tcd  s  (Oct  t  €,j)d¥<lu  (11-17) 

and 

Tu  * 


0 


(11-18) 
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If  Aq  is  constant  over  the  entrance  pupil,  then  £q.  (11-17}  is 

Sg 

where  3g  is  the  area  of  the  spatial  frequency  sons  associated 
with  the  spatial  frequency  Equation  ( 11-19)  is  a  useful 

tool  in  determining  the  transfer  function  for  a  diffraction 
limited  optical  system. 

The  transfer  function  which  may  be  attributed  to  the 
turbulence  alone  is  found  by  dividing  (11-16)  by  (11-19). 

Expressed  in  complex  form  the  transfer  function  for  the  turbulence 
is 

t-  Xc+iT, 

A.1  St 

The  concept  of  spatial  frequency  zones  makes  it  clear  that 
there  is  a  sharp  cutoff  to  the  spatial  frequencies  transmitted 
by  an  optical  system*  Since  the  largest  separation  possible  for 
a  pair  of  points  is  equal  to  the  maximum  dimension  of  the 
entrance  pupil,  D,  the  highest  spatial  frequency  passed  by  the 
optical  system  is 


f 


WAV 


This  is  the  Fourier  equivalent  of  the  usual  diffraction 
limit  expression  that  the  angular  distance  to  the  first  mill 
of  the  diffraction  pattern  is 


0C 


T> 


(11-19) 


(11-20) 


(H-21) 


(11-22) 


ft  p»rft<V«t  m . 


DtiMoywunT  o f  y<<Tof«lW  fatTots  fot  ftnuwllu  iyivaywiT 

diClOflifcYl . 


Wibn  ft  fv?ftTs  IV  on?  •  d’lvwtmitfYial  cas?  Oivxd  ssflibYift  Twoffs 


tV  lUio-dwAtmViwiAl  CUiY . 


$9fl\lVl  ft. 


'DwiOpVWfff  of  YQsToYQfttOYl  fatfOYS  £oY  0Y>(  •divWViaiOVial  T9YY>W*)I|| 


ivwomvff  diSfoYTiwi. 


lYnwioAiaTolv^  CoUiwin^  iWis  pq<|9  is  a  ^rApWttdl  YopioiQYiTrtTiovi  of 
tV Aitfotitm  ptortss.  TV  vmwtVw  of  divisions  was  <xvb(lV<xtili| 
CVftn  to  show  TU<  Aiff?Y9Yit«  in  1W  tffoCT  of  Odd  avid  W9n 


VwawVys  (&  dwisions- 
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6 


Original  image 


h 
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Distorted  image 


Graphical  representation  of  distortion  process 
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NoTtfriw 


F (q)  •*  Touv\cy  ttpttwtfMTtvvi  of  OY\<jW\  \vy\o^f. 

$(&)  '■  Vouy\oy  Y«^«wT<xTiw>  of  HisToYTQA  \vYift^«. 

yj.&)  :  FoiMlQY  YOfftttXffarfoVt  of  AVI  *lwi»XlT  of  H'lSfoYfad  inYHA^. 

Fl,  :  te  *tV»  ^’QYY'ft  of  V  (@^ 

?\i  :  Ir*TV»  +UYW  of  f(©) 
fe'iU  ttm  of 

6  •*  ravKj«  of  rc^ioYi  com'utaYod.  0  ^  9  <  2ir . 

A6  :  (\ft<^\ftY  St^AXftVlOK  of  «Ufrtt«VS  of  f(i*foYt«d  ivYlft^. 
v  A©  ■•  Av^uUy  di<t av\c«  of  0  dVvTovTion  cWvnoirtT  fyww  twtfwr 
of  d&wfiovt.  Y  *  s.  _t  o  ,  c  »  0,1,  ••• 

S  ■•  Kuwbw  of  diu&ttfvii  cf  dHtbtflon  (toftxxtbw  of  xUmxxrts  l««  f). 
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o,b  *  Ptatfav  co^’ic'iQYfli  of  9#*mA«K?A  bu  k, 


0  £  k  -  c  vi 


CV,U :  Foutwv  to$u\w(T$  o(  V\  iwAwA  bu  b, 


0  -  k  —  J^jlL  -  y\ 


VTO  Hutwboy  oV  'twptfT  po\v\tS  O^r  9Y\tWQ  iwML^t* 
m  :  !V\AqX  0^  C0t(4\c\QYft<S  ,  V\  ®  WH  . 


P:  P*-±-  *  Y* 


-  ,  »  5  i 

2S  Z 


cp  qj  -  2  p  ,.>r  *  0 


V  J-  >Y*  %  >° 


P  >  x*  O 
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F(e)  -  Z  k*0  (&w  tosh©  +1Sk^Y»feG')  CO 

UT:  F^C©)  5  ftfccoske  tB(,Swfe© 


1W*.  F*(©}  =  Xte  Tk60  (i') 

^(g")  *  ^Y.-l  ^(e  +  ^0) 

U)W«*.  ^(©  +Y&0^  '^\f.  p  t&W  CO*  k(©  *  V 

+  “Bk  s’w  te  ( ©  t  v  &G) 

Lets  ^ii  (©  tY&©)  *  p  E  Afe  COsfe  (QWbQ')  t 

+  Kk  s\y>  h(©  +V^0N) 

f  (0)  -  Lx  Z  k  (0  +  Y&  0) 

AssumVn^  twnvr»<|QvHf  o('  TU«  sew*  £  ( FouritY  SYYjes  ,  Hardij, 
l/Ul  WtllQtt  ,  I'HH,  p.31  ff.  ;  de  filAYltY  S»y‘w*  ftW<l  Ivnt^jU, 

GtmUiu,  T)ou»y ,  lYiVri  ed»fi»v»  t  IS'io,  p.TJO  ff.} . 
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f  (©)  =  Z*  Z  v  (©  +  v  A  6) 

L<t:  ffe  «  Zy  (©  +  V  A0)  (*) 

KSJ'Zkfk  ft) 

ExpOwAVft^'  Jb  (0^  s 

Sp{^k(cosh0  costa  A©  *s'»Y\k0  swta  A©')  + 

♦  Ub  (smk0  COstaA©  *  cost©  s'wfeY 

SUWWiVa^  (©  tYA©)  OV/QY  <  s  ”5.  ,  ,  SWHt 

5*w\‘X  s  “  Sin(*  X)  <xv\f\  sm  0  s  0  j  IVu  T«yw»s  (onTa'iyTtvi^ 

S\y>  tar&0  UhU  v&vusVt. 

Z  y  ^  (0  +  y  A©)  =  21  y  p  Cos  tat  A©  (  ft  b  Cos  I*  0  + 

+  'Sfe  sink©) 

Oy  fwwifsV.  fk?TfeI.,  p  cosfcvA0 
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SfcbsTiTiAttin £  \vi  W)  *?(©)*  21*  Xx  p  Cosfer  C±6 
It  is  now  w\A«tfT  tWf  isTta  fe-iVi  0?  'TV\«'Pi>ur»7Y 

S«t*S  f(0).So  ?rt>YYl(0  ftvui  lb<  IVlAvp«VlAftN< 

Ov  on\Ao<jofl*VfTi/|  0^  TUc  s'»y\  <Xv>^  cos  t«YW\S  •• 

<Xfc  2  I*  Xv  ptoskvAG 
kb  3  'Sfe  Zv  p  cos  It*  A© 


Tbtv\  •• 


ftfe  =  _ ®H - - 

21  v  p  Cos  bvA0 

s  bt> _ _ 

21 1  p  cos  b  <A© 

TV»«  lad  wpwdom  cu*  be  sswiplifWd  for  covnpUTafion  bi| 
folding  4Vie  ^utnwi atic>v> ,  since  Cos  (x)  -  COS (-  i) 
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prob\m  ot'oAA  t\v\A  oo 9^  nuwibwoV  ivtfwols  Utatefyi  cw 
d  4to  cwA'tnw  y  ^  0  ftv\A  subit’iTuriow  cj 

f  oy  p .  (  S««  notation  list) . 


T>w\opvn  rrTf  of  YoVToYOtf\ovi  fonTovs  fey  tuo-diw?Y>ViOMl  flyrtgoYotl^ 


ww&v'tarif  diflbYTiovi- 


LvY\wcA\ttT«i^  foUouhn^  IVus  sIa«$T  (xyc  <jYapvt\((Kl  YQpveswfldtiofls 
of  iV*c  t«*i|)0Y<kU^  iway'icmT  d'trtotfiw  prows.  lv>iVm« 
W\Yf«s«<rtdT\ovtt  fy(U  s^uoor*  Y«pYos<mfs  (\  po\y»T  valu*  of  1W 


!wa^c  'TVi«  v^iayoIooy  of  dii/isiow  U)as  OrtbilYAril 

fcbown  to  s\aou)  ddfwvu?  \i\  tffcfl*  of  OdA  Gvwl  w<n 


1 


v^bm  bf  dwisiow. 


SIO  Ref.  63-10 


-91- 


i 


-ZA£-SAfl  -Iftfi-Afl  4fi.  JAfi  3Afi  ZAfl 
2  2222222 

Distorted  Image 


Effect  of  folded!  summation 


Z 


r,8 


.  49 


r  -  -  m  i  o 

s  *  -nio 
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H<JT«T\OV\ ; 


C,T),l=  ,F  *•  Fouwr  <ot(nc‘t«<lTs  u)Yft\  /9jp?(T  To  "S(©,^) 
C,d,«,?:  Fow/ioy  to<W»it*ATs  ujVtVi  y? sport  To  €  (©/£) 


:  av»  rtoroodT  ok  "Wm  distorts  i 


« . 


X  ••  IWiT  Ob  X  (h  wmboi  o(*  tows  of  To<ayV?y 


-J  :  1\^aVT  ob  a  of  v\vA«\b>’v  of  f»Y<ns  of  Fowriot  for  Vs. 
c  s  ivwUx  ov>  X  ot-  f*Y  *b*  ,  L »  0, 1 ,  •••  ,  X 


J  '■  iM>'y  Ob  ^  ol*  ,  j  -  0,1,  •••  ,  *0 
1*1:  biAwb**  of  d’tvMtObs  of  ciiiToTtifcw  owK.'lWf*  aw  Mti 
N*.  v*w«ibot  of  d’vjiswbs  of  disTotliob  6b  j  .X'o'tbTs  ax*  N+l 
Ww  vndox  Dm  a  OV  ^o\bfs  of  ditto  flion  ,  w  -  -  K 

0  :  iwAu  Ob  of  poibfs  <jf  distortion  ,  b  =  -  N_  K 
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p  *.  UiO^VlTrt^foflOY  ;  ps|JY#i:il.>S^±^. 

*  -t  i  v  *  ±  £L  ,  s  *  ±  4 

=  J-  ,  Y  *  £  $_  ,S  =  ill 

»J-  ,  X  *  ±-4-  ,  s'*  il 

c| »  w«\(jVvT^  ^  tatfoy  *.  c|  =  *f;  y  f*  0 ,  &  i  s  ^  0)Jl 

-  1  i  Y  *  0>*L  .  S*  O,!! 

2  2 

-2  i  X  *0,  JUL  ;S*  OjH. 

2  1  2 

-  •  >  <  *  0,  &  .  s*  0,SL 

2  J  2 

ti  VAUH\pli9Y  ft©  :  yA©  is  AV^QvaUy  dlSlAVWC  Ow  %  0 f  ?|«W»AT 

Ww  CtVtfoY  o(*  ^SfOYtiOVl.  Y  -  *  o 

2 

<5-.  vNttKi^ri«y  Of  Ad?  :  sA^  U  Oin^uUr  disfoo**  on<|  o£  el»fo?rfT 

(yoyo  c?Wh<  of  dfetlrtfuM.  s  *  JL  -  n  >  0 

z 
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*.  ObscilSfc  o(  field  . 

^  ••  ovdivtaTe  of  iw\a<[<  field¬ 
's  (0>$  :  Fouy’hy  Yopm**Tafiov>  of  undiflctfftd  iwa^. 

2  (©  •  TouyIqy  YQpTtstnTaTibn  of  dht6rT<?<i  i w\  . 

@  :  M^Uy  W1«&UY«  ovi  ^  o('  field 
<t> :  avj^Uv  vm9«suy<  oy>  of  field 

G\vj»w 

.<£)  ‘^cm  COSL  0  (OS  j^)  /DijCOSt0si«j4>  t 

+  £i.j  Si*u0  (OSj<t>  +  F;j  sin t6  S»V1  j 4>) 

L«T*  H^(g,4^  *  Ci]  cost©  ccsj<t>  t'Uij  COSL0si^j4)4 

+  £  i]  Si« i  0  Cos  t  T L  j  Sim  0  si  n  (|) 

--  2L;  Xi  £ij  (e.4>) 
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/  A  11  ft 

z(0  »<WS  Z.*t  <^(©+Y£s©  ,<j?  4  S  k4>) 

Ulortt*.  ^(0tYA©,<(>  * 

"  ZT  t  Z  i  plCy  '*■  (©  4  Y £^6l)  Cos  j  (<t>  4  sM)  + 

+  (DSt(©4  Y&6i;S‘i/^(4>  +  4 

4  s’»rtl(©  +  yCs6;  Cos^(4?  + + 

4  KjSWi^©  t  YCi©J»Si*  j  (<p +SA<£) 

l«t  ^.(et Y^e,<5 1 *  pfcj cou(e+<^e) cosj (<t> ts&$)  + 

•4'Dij  (Qu(Q+XL&)%\y)j($>  +  sL$)  + 

*  £cj  Stfll  (01  Y  A6) cos] (45  i sk$>) 4 

♦-^L]  $'»*i(04  Y&ejsifl]($4sC^)  (2) 


z  '  ZvZ*ZiZ.i  (etxae.ots^ 
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ftssuwivu^  lo*w<!Y<^ftce  of  s?ytos  as \vi  Action  ft 


*  *  S-ili  (l, 

L*fli  “H  =  2Ly  X*  ^ ^ j 


'Dft’n:  -  £;  Z\ 


W 

(H) 


£xpaot\t  (z)  ••  +  +  = 

’  pfclj  (cos  t0  CO.S  lvAG  -  sVni.0  S>n  cx60J* 


•  (cosj<t>  coi]s A4>  -  swj<£  sivijsk^  t 

**£>;•  ((DStOCOSLXft©  -S\nl0sinM©V 

•  (s\*j4>  cos^s&4>  +  cos j4>  s'wjsM^  t 

+  £ij(sinl0  coSuA0  t  COSt0Si'OtYft0> 

•  fcosj^  cos  *"  sin  -t 
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+  Vtj  (s'mi©  cosixk©  t  COsLSsmW  A©} » 

•  (sw»j4>  cosjs&4>  +  cosj<i>su>jsA<t>) 
taboj?,  ou*y  Y  =  -jj.  ,14.  , 
s  s  -  K.  j  il  t  cotftlM sivi  u&O  owA 
s\rt}sA$  o)»u  vaa’isW  s\*t?  s\v3  *  *  -sin(-*)  avu\ 
s \n  0  -  0- 

ZxXs  (©  4Y&G,<t>  )  =  ZvX*  p  COSIX  A0  COS  js&$# 

•  [cij  Cos  i  0  cosj  4>  t  x»l]  Cos  t  9  SIV»]  <)>  + 

♦  ELjSinLd  cosj 4^  f  Tij  S’miGSm^J 

Fv»w  (i)  OwA  (i) t  =  £-  XxXs  p  cosivAe  cosjsA^, 
Gwd  subfcTiTtftM^  In  Qi)  • 

*(M>  ZtXi  p  cosivAG  cos  js&$ 
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ZttSAM  WidwT  TUcTi  2.1'^  is  The  oHv>» 

Toutin'  Y«f<ts<flfifl^  l(0,4O-  Twwilto  \v\dyfldotu? 
OY  0^1\A«^OV\ffilifi|  ol  tOWpBwAf  f »tfOS  : 

^  Li  '  _ ftii. _ _ 

Iv  li  pCOSW  COs]sA$ 


<v- 


Zv  "Zs  p  COslvkO  GOS^S &4> 


avid  l»h<?uAs*  fov  rD,&.,a*dT*  ujTtTVi  Y?sj><?<T  To  d.savidV' 


oU  t  j  . 


fts  i*  2<?duvi  ft,  tW  doviovrt‘»rtdToY  abov?  Cav>  b?S*itApl»ifi?d 
(ox  (.OM^uUfiovi  bo|  foldiv*^  iVi«  SummaTioms,  siv»(f 
Cos  X  •  CO  S  3C)  .  TVic  pfolol'WI  ck  odd  av»4  «v*n  iwfwuato 


is  tab**  cav«  b«j  ito  covidifUiAS  V  ^  0 ,  s  ^  O  avid  b 


H 
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subtfiTwtiDv*  oP  Pov  p.  Tfo<  IVnyi  bnows: 

Jl  11 

2Lw*o  2Ln*o  c|  cos  U  A©  Co*jsA<t> 

U)W«.  X  *  Jl_  -  v*  >  0 

s  ’  ji  .  n  ^  o 
2 

W  ii'iOkWt  H  aAWy  IftiUwdttV) 
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‘T&£Tc>r<rriffYi  fct  Yaviftyfl  Hitlwriow  ■ 

ft  tW  Ow-divYiovttiew)  Cos?  ftvwl  Srtfl»W*&  tftftfs  IVk 


tuJo-dWv\t\*W  OHM. 

Snt>fA  ft  . 

■foSwort’iw  J>V  ov\«  *  dWv\iioM\  iwa<y  U1\T\a  TtvnpOY^lvj  Wiftnl  (torJUn. 
G\v)«V»  :  ft  V\IAV*\W  o(  d't{fotf»A  lvY\tt(|«i  (tfdwd  j  *  I,  Z,  4"  ,  «atV» 
vtprttwfbA  Ok  fouvw  soy'ws  gJiTVY  cotff »ci«Y>Ts  At  >j  ftv»d  1i;,j 
ividowA  w(TV\  v«s^(T  to  fvtauwuij  fe>^  i  *  l,!,’*’  . 

To  9\iu\\  Tit*  COtffu*»«>rtTs  ,  Aft  GVUi  'Sft 

(Wrfiwv- 

X  l  1.  ^  \ 

(I)  A*  8  VY»a*X^  (ftouj  ) 
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(z)  ft  ft,  *  £  \  ftg.L 
A  2l  ±  ) 

(Z)  (&*)'  ^  ( £  i 

(4)  (■&»')  5  ( z-  "B. 

SoluTiovt : 


ft  ft  *  Zj  A*^  j  ?Y(W»  (t} 

£j 


Su\&£fauf»v\^  \v\  (  0  ; 

'Bo,  Z.i  ft  <M_ 

riTSft.j 


t  13a  '  VW0Cj(A*,i  +^A.j  ^ 


*Ga  f  £4  ft<M  +  1 


WfcXj  (A4,;+^a,n 


^4 


1  ~j  ’4 

L  gVj^i  H 

2il^i 


ft 


T^a  £j  ft*»L 


Z  i  ISft.j 


•  s^v)  (Z  j  A 
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^cHgjQ  « 

&  TUtt-dwnmio'vml  iwa^c  uiCCVa  YaYiarif 

dstottvovN . 

G'\o«v>;  C\  mn\ow  o f  diflwf«A  lw\a^  iwdwd  Wv <  1,1,  ••*  ,  oad> 

Vt^Yfod  Ok  Youviiv  S<J<i*s: 

fw  (ek,4>0  =  Z;  Zj (citi«c«s iei.  cosj^!  + 

*  'Dt.jy*  cosce^ s\y)  j<t> x  ^£^wsmtekcosj41 t 

+  F^w  sine^si/i]^) 

TW  <^\/9y\  ^u<0Uy)  com  b«  Y<?ujv’r1t<?v>  as*.  . 

^w*(©|i  i^l)  *  -L-  c  X j  Cos(t0|,  t  j  4^x)* 

+  Fiiw)  Cos  (c0it  “j  <0  ♦ 

♦  foijw.  +  Eij„)  «'i«  8  k  t  j<f>i)  t 
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Sivwc  IV  vesTottfrion  pfoc«js  is  IV  fov  tocli 

t  owA  j  viffTCCriovk  is  dvopp<?A  »v\Twc  following 
Tof»\ftA:  TV  /QiTor<t4  Cotf£i<i?n1s 
Oovk\iTisv»$  •• 

(100)  (C-fV  +  (W)'  =  max„  [(<?„- FjVou +  £«)*]  *  (*!) 

(101)  (OTV+fatey*  [  (o  F„y  t  C-D„+e„j]  •-  (**) 

(101)  C  - F_  *  Z»  (cw- F»)  =  (*3) 

T  +  £  Z  w*  (T*>  +  £*0 

Ool)  C*  »X  «  (CmtrJ  ,  (#4) 

T>+£  Zw)  (-Pw> 

C»OHj  SV(C-rt  *  ^  (cfyy,  'Tw»)] 

Oofj  5^  (c tr)  *  S<^w  [x  +7V)  ] 

(\°0  S^n  I'D  t£  )  *  [  Z.*v»  (l)w  +  CE.  r*)J 

(107)  =  s%»JjCvr>  (-Vnt£*,)] 
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<j»..l)  Ftwvv  ( loi)  ■■  (C-FV*  (t)  +£■)’•  (♦S')' 

(102.1)  Tt.««  (lol)  ••  (i  t  Fy  «  (-  »  *e]'  ■  (tf  4)' 

SuWl'itiTioO  ( I0Z .i)  in  ( too)  aw  A  (lovi)  in  (lOt)i 
Ooe)  (d-hOM**}’ +  (utey  - 
(ioiJ  (-T>tE.y(**»V'  +  (-X>it')x-  (*£) 

Fxwn  (too)  anA  (lOt»): 

.^n[l«  fa™ +£»')]  •  (*£) 


(uo)  Co  +  E'i  « 

(#i) 

Fyow  (ioy)  ('07)  • 

(u»)  (-d+eV 

(#  0 

.(*4)ltl  . 

•  Y1 


[2«.(-D* +£«)]*(*  8) 


lloi.l!)  Ftww  (to*)  :  (D+eV  -  — L-  (C  -  r)^ 

(*t>y 


-A2- 


(kh)  F«m  («oV)  :  ( ""0+  £)'  *  JL  (  C+f'j 

(#4)' 
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SttbsfiTifiW^  (loi  ‘"l)  \v\  (lOO)  : 

(JI2J  (0-F')'4  J-  (<!-F)a*  (tfi) 
(|03.'<i)  1\A  (loO  • 

(M3j  ((JtFj'i-  J  _  ((Ttn1-  (*Z) 

(*«0 


Fy6*>  ( IU)  a/wi  (l04)  ! 


O'*0  (<T-F)  s 


l'1Lsy«[2;«('C'm-Tw)]  *  (*«) 


Frfcw  (  U3)  (xf\A  (lOt>): 


(Utf)  (C+T)=  — —  v  .  t  Fwj^J  •  (*7) 

L  (#5)'  +  1  i  J 


Ad<\w»<i  (ll**)  Uvwl  ( ii(>)  : 

C  =  (*£}  *  0*  ~Q 
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Or .  C  «  _L 


VYlWy*»t(^W»' Fwi)  ‘ 

fr  (*Dv*  t  £*n)  1  ] 

Iw^w't  £m) 

>•  1 

W\Wv* 


[(Cyy\  t Ttft")  +  ('Dm  +  £w)  J)  * 


Zw.f-DmtEw) 

^ T w”) 


*  \ 


.  S^Ive(C»etFtn)] 


d  -  (*  6) :  (* 

7. 


Oy:  De  -j? 


+  (Dw  *  £*»)  ] 

|Vk 

\Zm(  Cm  'Fvw! 

+  1 

* 

L  .  Z.w»(T>w»  +  £»n). 

1 


V 


VttOflvn  [(CmtTwTt(-Dn,t£wT]'| 
'l»fatFw)  ]*  -.  I 
.  Z.  *e(‘T)w+  .  J 
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A44in^(ll6)  avid  (in)  : 

£  -  (*£)  +  (*8^ 

2. 


Or  E  -  _L_  f  Vvxaxtm  [fcw-T^  ^(x>\wt£\w)  11*  _ 

I  ZwC^w-Tvw)]2  4-  |  |  .S^ik^fDwi+£w))  + 

t  .  Z  m(Dmt£***y  J 


+■  If  fftkx  w  Fw^  +  (~T)»>>t  £v»^  ]  )  ^  -i 

fx. w.  (Orv*  *  4-  I  .  S<jVi[l*,(''Dw*+  E  y*t)J 

LX  vn  ( Etw^  j 
SubTrutTtv^  ( II 4)  fww  (ll^ ; 

r  =  -  (ftg)  -t-  (»7) 

2. 
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Ov  F  *  «  J_ 

t  (t)v. + £~V  ]  1 

l 

|  (*D>w  *r*\ 

1  1*  1 

\  1  ?  wi  (C\*»  '  T v+)  , 

4 

J  VAAXvw 

’(C-t7-T+K*E«)’J  )SI 

|  rZ.J-'Dw*t6*} 

'•'  i 

l  „Z wi( Cwi t 

.  s^n[  2  ««((?«•*?<*)]  + 


W  illiftm  Hadbu  IftitaMton 


